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^ . Abstract: We demonstrate the effective action scheme for the leading parton 

interactions and discuss the symmetry properties. The interaction kernels are 

CN I particular cases of conformal symmetric two-particle kernels. There is a direct 

' relation to the conformal symmetric rational solution of the Yang-Baxter equa- 

'sj" ■ tion. 
(N 

(N 
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-5 ; 1 Introduction 

Ph. 

Ph. During three decades high energy processes related to the Bjorken limit of scat- 

r-| ! tering amplitudes are playing a major role in the investigation of the hadronic 

structure and the interaction of hadronic constituents. Quite a number of theo- 
retical concepts and methods have been developed resulting by now in a standard 
treatment described in textbooks [1] and reviews [2]. 

In the last years increasing attention is devoted to topics going beyond the 
standard deep inelastic structure functions, among them the long standing prob- 
lem of the explicite treatment of higher twist contributions and their scale depen- 
dence [3], the generalization to non- forward kinematics [4], interpolating between 
the scale dependence of structure functions (DGLAP evolution [5]) and the scale 
dependence of light cone meson wave functions (ERBL evolution [6] ) as two 
limiting cases. The existence of this interpolation has been pointed out early in 
[7]. 

The small x behaviour of structure functions and related questions have drawn 
the attention also to the relations between the Bjorken and the Regge limits. The 
evolution in the latter asymptotics, as far as it proceeds in the perturbative region, 
is represented by the BFKL equation [8, 9, 10]. 

In view of these topics a discussion of the Bjorken asymptotics invoking con- 
cepts not aligned to the standard approach may be of interest. We advocate 
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an effective action approach [11] in analogy to the one developed for the Regge 
asymptotics [12]. 

The light ray position as the essential variable of operators in the Bjorken 
limit appeared explicitely first in the approaches by Geyer, Robaschik et al. [7] 
and by Balitsky and Braun [14]. 

The parton interaction operators can be obtained as external field effective 
vertices, where two type of external fields are introduces describing the asymp- 
totic interaction with the currents of high virtuality and with the hadrons. This 
calculation can be done analyzing the space-time field configurations only without 
transforming to momentum space. 

Prom the very beginning conformal symmetry was underlying the ideas about 
Bjorken scaling and its violation. It turned into a tool for finding multiplicatively 
renormalized operators [6, 15] and for relating the forward to the non-forward 
evolutions [13, 16]. Combining conformal and supersymmetry leads to interesting 
relations between different evolution kernels [13]. 

Integrabihty of the effective interaction in high energy QCD amphtudes has 
been discovered by Lipatov [17] in the Regge limit for the case of multiple ex- 
change of reggeized gluons and a similar symmetry property was expected in the 
Bjorken limit. The first examples of higher twist evolutions tractable by integra- 
bihty have been studied by Braun, Korchemsky et al. [18] and by Behtsky [19]. 
Relying on well known methods [20] solutions of the Yang-Baxter equation with 
conformal (sZ(2) and superconformal (sZ(2|l)) symmetries have been ontained 
[21] in a form suitable for the Bjorken asmptotics. Some results will be reviewed 
in sections 3 and 4. 

2 Parton interactions in QCD 

The standard process is electron - proton scattering ep — > eX at high energy \/s 
with a large momentum transfer q, —q^ — xs, x — 0{1), from the electron to the 
hadron. The hadronic part of the process can be written in terms of the forward 
imaginary part of the virtual photon - proton amplitude. 

The non-forward deeply virtual Compton scattering off the proton , 7*(gi) + 
p{pi) 7* (92) +p{P2), is an appropriate generalization (Fig. la): 

qi/2 = q' - xy2p', q'^ = 0,p'^ = 0,p' ^ pi, s = 2p'q', (2.1) 

-qi/2 = Ql/2 = xi/2 s. 

The Bjorken limit corresponds to s — > 00 with the Bjorken variables Xi and/or 
X2 being of order 1. 

Let us represent the light cone factorization graphically as in Fig. lb where 
a sum over t-channel parton intermediate states is implied. The approach to the 
asymptotics can be visualized as a process of iterating this factorization. The 
next step is shown in Fig. Ic. 
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Figure 1: Large-scale factorization of the amplitude. A summation over t-channel 
intermediate states is implied. 

The Green function H (not necessarily connected) represents the parton in- 
teraction. The effective action is to describe just this parton interaction. 

In the lowest order of perturbative expansion there are only pair interactions. 
The partons can be identified as modes of the underlying gluon A and quark 
/, / fields in the gauge q'^A^ = after integrating over the redundant field 
components p'^A^ and p'^j^ ip. 

We represent 4-vectors by their light cone components x± and a complex 
number involving the transverse components x_i_ — Xi + 1x2- In the case of the 
gradient vector we change the notation in such a way to have = = 

l.d^x = d*^x* = 1. In particular the complex valued field A represents the 
transverse components of the gauge field. We choose the frame where the light- 



like vector q' has the only non-vanishish component q'_ — ys/2 and p' the only 

non- vanishing component p'_^ = ^J s/2. 

We decompose the Dirac fields into light cone components, 

il) = ip_+ 7_^+ = 7+z/'_ = 
s 

Ua,b, a, 6 = ±, is a basis of Majorana spinors, 

l+U-,b = 7-1i+,6 = 0, ^Ua - = l*Ua,+ = 0. (2.3) 
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The gauge group structure of the action has been written by using brackets 
combining two fields into the colour states of the adjoint (a) and of the two 
fundametal {a and *a) representations: 

(A^TMs) = -iP^^Af A\, {jlf^h) = t%!T!L 

ir^A) = t'^pJ*^A^, iA*t*-f) = t'^^A*'/^ (2.4) 

We shall restrict the detailed discussion to pure gluo dynamics, presenting the 
final results with fermions included. 

The action in light-cone form can be recovered from the kinetic term 

= Jd^x- 2>l*"9+a_A" - d^dA*''d-^d^*dA'' (2.5) 

by extending the transverse derivatives in the second form of eq. (2.5), relying 
on the residual gauge symmetry, 

d^A" {V^AY = d^A" + ^(A*TM), (2.6) 

The result can be written as 

S = J d'^x{-2A^*{x){d+d- - d^dDA'^ix) 

+l{diV^23[diM^i)iA*ix2)T-A*{xs))+c.c.) 

+^ Vii'M'^iA''i^i)T'dA{xr){dA*{x2)T'^A{x2'))U=,'=. }. (2.7) 
The elimination of redundant field components has lead to non-local vertices, 

n*23 = - ^3) + cycL], Vn',22' = (^1 + dr)-'. (2.8) 

Here and in the following we omit the space index + on derivatives, i.e. derivative 
operators not carrying subscripts — , _L are to be read as 9+. Integer number 
subscripts refer to the space point on which the derivative acts. The definition 
of the inverse is to be specified. 

3 Space-time picture 

The block H of the virtual Compton amplitude describes the interaction between 
sources located in the vicinity of the light ray x± = 0,x_ = 0,x^ = z E 
and other sources the distribution of which is almost constant in the variables 
x'j_,x'_^ depending essentially only on the coordinate along the light-ray x'_ = z'. 
The small width of the first distribution near the light cone is characterized by 
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the short distance scale A ~ and the variation of the second distribution 
in directions away from the hght ray is determined by the large distance scale 
5 ~ m~^. 

Consider now the QCD functional integral with such sources or the related 
vertex functional with corresponding external gluon and quark fields, we di- 
vide the fields of quarks and gluons into two types of external fields A^^"* and a 
quantum fiuctuation Aq, 

A^ A^+^ ^ Aq + A^-\ (3.1) 

has the support near the light cone and has to be substituted by the following 
expression in terms of (regularized) delta functions, 

A^+\x) = A^{z)5{x^)5^''\xA_) + e»(A). (3.2) 

The other external field has to be substituted as 

A^-\x') = >li(^') const + C>(m). (3.3) 

The vertex functional or the external field effective action is now obtained by 
doing the integration over the quantum fiuctuations Aq. Consider in particular 
the resulting vertex involving two A^"*"^ and two A^~^ type fields on the tree level. 
It has the form 

j d'^xid^X2d'^xi>d^x>A^+^ {xi)A^+^ {x2)G{xir)G{x22') 
[Vli'G'(xi,20l^22' + S^^\xr2')Vir,22']A^-\xv)A^-\x2>) (3.4) 

Contributions from disconnected graphs have to be added; they result in the 
cancellation of IR divergencies. 

Vw, V22' are simply related to the triple vertex in (2.7) depending on the case 
considered and G{x) stands for the quark or gluon propagator. 

For example, the triple vertex contributes to the tree graph Fig. 2a as 

and reduces in application to A^^\xi), A^~\xi') to 

Vir = -tg I d'x,x'xv-^^^'-^Al{AlT'^Ar)U,=,,, (3.6) 



di + dy di 

We substitute the particular asymptotic form of the external fields (3.2, 3.3) 
and obtain 

Q'^ f [~ ~ 1 

cln— dzidz2dzrdz2'AiA2ArA2> V1VV22' J a + Viv 02' Jq (3.7) 

5 
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Figure 2: Triple and quartic vertex graphs. 

We abbreviate the residual dependence of the parton fields on the light ray po- 
sitions by indices 1,2 for Ai{zi), A^^\z2) and by indices l',2' for A^^^'{zi>), 
Ai ^\z2'). The integration over the transverse and + components of xi,X2 is 
done due to the delta functions and the integrals over the transverse and + com- 
ponents of xr,X2' are represented by J a, Jo, 



cln — -Ja 

Q 



dxi'+dx2'+cfxii±cfx2'±d^di* (x 



2 \-l/™2 



11' 



"22' 



X' 



v-1 



1'2' 



cln — - Jo— dxi>+dx2'+d^xi>±d^X2'±didi*{xli,)~^{xl2>)~^5^'^\xir2') (3.8) 

777- •/ 

The integrals are regularized by taking into account the smearing of the near 
hght cone distribution by A and the large scale cutoff S for the other distribution. 
We shall do the integration in the logarithmic approximation. For this we have 
substituted in (3.7) already the effective triple vertices omitting terms which do 
not result in logarithmic integrals. 



The integral over the transverse coordinates is logarithmic in the region Q~ 



< 



\x\ P < m ^. 



1 J 
cln— Jo 



dx'_^_d'^x'j_S{zi' 



[-x'+zii' - + ie){-x\z22' - is) 
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/ TTI? / da{azn' + 1 - iey^{az22' + 1 - ie)~^S{zv2') (3.9) 

JQ~'^ \Xj_\ J—oo 



The result can be represented in the convenient form Jq = Jq^\ where 

Jo'\zi, Z2\ zi', Z2') ^ da Xo'\a)S(zir - azu) 5{z22' + (1 - a)^i2), (3.10) 

JO 

The case of Ja is a bit more involved and can be done as follows. 

cln — -Ja = / dx+iix+2'd^x±iid']_2'{xi' ■ 2;^)(— x+i/^^n/ — \x^f\^ + ie) ^ 

{-X^2'Z22' - \X2>\^ + i£)'''^{x^V2'ZV2' " kr2'l^ + '^^)~^ (3-11) 

We use the Fock-Schwinger method to do the integrals over the x^ variables and 
obtain 

daida2da^5(^ai - l)aia25{aiZiv - a^zy2')5{a2Z22' + Oi^zv2') 

j d^x±_vd']_2'{^v ■ ^2') {(y.i\xy\^ -\- a2\x2>\^ -\- a2\x'\;,2i\^)~^ . (3.12) 

The logarithmic contribution of the transverse integral results in 
j daida2da^5(y^^ai — 1) ^'^^ 



(q;2 + a^){aia2 + a2a^ + a^Oiiy 
5(aiZn> - azzv2') 5{a2Z22' + o-2,zv2') / TTTF- (3.13) 

The convenient representation of the result is Ja = Jiii, where 

r(ni + 712 + rig - 1) , , , r/Y- 

r(ni)r(ri2)r(n3) Jo ^ 
<-^q;^^-V^3-' 5{zn' - aizi2)5{z22' + 0:2^12). (3.14) 

The general definition of the a integrals appearing here is 

— [axi + 1 - ie] "1 [q;x2 + 1 - ie] 

-00 ZTTZ 

...[ax^v + l -ie]~"'^. (3.15) 
rii are positive integers. The Feynman parameter representation is 

r(ni)r(n2)...r(njv) i 

5(E«i^^) (3-16) 



In the latter integral the varibles ai range from to 1 and V^^^a — da\...daN5{{l— 

By decomposition into simple fractions one can reduce the number of factors 
in the denominator of (3.15). In this way one derives relations like 

Jl2(xi,X2) = Jii{Xi,X2), 

Xi - X2 

Jui{xi,X2,X3) = — Jii(a;i,X2) Jll{xi,X3), 

X23 X23 

Jll2(xi,X2,X3) = ^^^^ Jll(xi,X2) + ( — + — ) Jni(xi,X2,X3), 
X13X23 \X23 X13J 

J22l{Xi,X2,X3) = —Jui{Xi,X2,X3) 

-^(|^ + f^)^"fe.-^)- (3.17) 

2^12 \X23 X13/ 

Wc consider now in some detail the effective interaction of gluonic partons. The 
vertex in the external field effective action A*^~^^ (xi)A*^^\x2)A^~\xii)A^~\x2') 
referring to the case of parallel hehcities (chiral odd channel) gets a contribution 
from the original quartic vertex (2.7) and one from contracting two triple vertices 
by integrating over the quantum fluctuation Aq. The triple vertices in the Bjorken 
limit corresponding to Fig. 2a, b are 



^2 ^*J_ ^2^J_ 

(3.18) 

Writing down the extenal field vertex as in (3.4), substituting the asymptotic 
form of the external fields and doing the interations besides of the ones in the 
light-ray coordinates we obtain as the particular case of (3.7) 



/ 



dzidz2dzi/dz2' 



djdl,+dldl did2' + dvd2 

(91 + 910291^2 '''^ (9i + 9i02 ° 



{AIT" Ay) {AIT''A2>) (3.19) 



The first term in the brackets is obtained by contracting the triple vertices by 
substituting the points 1, 1' in the second and 2,2' in the first vertex in (3.18) 
and vice versa. The second terms emerges from the quartic vertex (2.7) for the 
case that A\^ A2 are of A*^"*"^ type and the remaining two field of the type. 
The second relation in (3.17) Fourier transformed to light ray variables implies 

where 

jg^/ da^^^^dizn'-az,2)Siz22'), (3.21) 
Jo a 



8 



X^^ = -(1 - «)^ xS^^ = -(1 - «), x?^ = «(1 - a). 

We use this relation and the one obtained by interchanging 1,1' 2,2'. 

The remaining divergence signalled by ^ in the latter integral is cancelled after 
including the disconnected contribution of order to the self energy in the two 
propagators connecting A\ with Ay and A2 with A21 separately, 

2 j dzidz2dzvdz2' Wg 5{ziv)5{^22') {AIT'' Ay) (A^T^As/), (3.22) 



'^-Mio^t ^^^^ 



= --(- - wf = (3.23) 

Adding (3.22) to (3.19) the singular part of Wg results in replacing in j'fy (3.21) 
a by adopting the conventional "plus prescription". In the following we 

understand the factor ^ in Jiy with this prescription. 
We obtain 

j dzidz2dzydz2'[4f + 5^') + (1, 1' ^ 2, 2')] 
{d^^AlT'^dAr) {d2^A;T''dA2'), where 5^2) = s^^u') S{z22') (3.24) 

We have calculated in the space-time effective action approach the vertex of 
parallel helicity gluon interaction to order g"^ in the logarithmic approximation. 

In the case of antiparallel helicity gluons the calculation of the effective two- 
parton interaction results instead of (3.19 ) in 

\{di + dvYdid2 {di + dvY 

-^llw^ ^0 iAlT^A2) (AyT^A;,) 

+^^^^1^-^111 (AlT^Al,) (A2T^A2>) (3.25) 
O1O2 

Here the self-energy contributions are not included yet. We separate in the first 
bracket a term equal to the result for the parallel helicity case (3.19). Adding now 
the self-energy contribution (3.22) results in the replacement of the first bracket 

by 

2dyd2'[Jn' + «;f - {dy - di){d2' - 82) Jin + ^i^sJo, (3.26) 

where now each term represents a regular operator. 

The last two relations given in (3.17) for the J integrals Fourier transformed 
to light ray variables can be used now to substitute {dy — di){d2' — 82) Jm and 
{dy + di)"^ Jill by J112 and J221 plus remainders involving Jq. This results in 

(2[Jn' + 5(2)] + j^^^ _ 2j^^2) {d^^AlT'^dAy) {d^^A2T''dAl,) 



+J221 (d^'AlT'^dAl,) {d^^A2T''dA2') 
+jQ^^Q^Jo {-{d^d2' + d2dr){AlT'^Al,) (^TM^O 
+{didv + d2d2')[{4T''Ar) (AsTM^,) - (AtTMs) (A^TM;,)]}. (3.27) 



Due to the commutation relation of the generators we have for the gauge 
group brackets (2.4) 



{AIT" Ay) (AaTM;,) - (^^^"^2) (ArT^Al,) = (^*rM;,) {A2T''Ar). (3.28) 



Thus the last term in the brackets multiplying Jq in (3.27) is equal to the first 
one up to the sign and the exchange of 1' and 2'. We remember that the above 
expressions are understood as integrated over the light ray positions 1, 2, 1', 2'. 
Therefore the exchange of 1' and 2' is just a substitution of integration variables. 
We conclude that the contribution involving Jq cancels and arrive at the result 
(3.31) up to normalization. 

The resulting effective action involves two types of fields A'^{z) hving on the 
light ray. The (non-local) vertices in this action are converted into Hamiltonians 
by considering A~ as annihilation and A^ as creation operators of partons. The 
Bjorken asymptotics is calculated by studying the evolution generated by those 
Hamiltonians in the (Euklidian) time variable ^, 



We list the operators of the leading parton interaction in the QCD Bjorken limit, 
the momentum kernels of which are well known, in the light-ray and Feynman 
parameter form. We restrict ourselves to the ones involving gluons (A) and/or 
only one flavour and one chirality of quarks (/). We use the abbreviations for 
the kernels (3.10, 3.14, 3.21) and also — Z12J. The position arguments of 
the fields will be abbreviated as indices 1, 1', 2, 2'. We suppress the label (±) 
distinguishing creation and annihilation operators, since the fields at points 1, 2 
act always as creation and the fields at the points 1', 2' always as annihilation 
operators. The integration over the positions, summation over colour indices and 
operator normal ordering is implied, 
parallel helicity interactions 




(3.29) 



{4[jif, + (5(2)](a-^Atr»MiO - 2t[4i} + wfS^^^]{d-^f*t''fr)} 

[-2(9-U;™A20 +^(9-V2^"/2')] 
-Ai4^,{d-'f*rdAr) {d-^A*t*''f 2') +C.C. (3.30) 



anti-parallel gluon interactions 



{8[jif, + 5(2)] + 4J221 - 8 Jn2} 
{d-^AlT''dAr){d-^A2T''dAl,) 



10 



+4J22i(9-UtTOA*,) (9-1^2^^20 + c.c. (3.31) 
anti-parallel helicity quark interactions (one chirality, one flavour) 

{_2[j(f) + ^(0)5(2)] + Jni}{d-'f:efr){d-'hef^,) 

+2ji'\d-' h){d-' hr f;,) + c.c. (3.32) 

anti parallel helicity quark - gluon interactions 

-Ai[J^, +wf5^''\d-^AlT''dAv) (9-V2i"/2*) 
+ (5(2)](9-Vi*r/iO(9-iA2T"M;0 
-AiJ2ii[{d-^nt'-dA\,) {d-^A^r^h,) - id-'AlT^dAy){d-'f2t''f;,)] 

-8iJiii{d-^AlT''dAr){d-^f2t''f;,) + c.c. (3.33) 

annihilation-type interactions 

-4iJ{u {d-^Airf*,){d-^A2t*''f2') + 2i5(2)- (^-^^^^"^^(^"Mat^/a* ) 

+12^j(^)- {d''AlT-d-'A2){fyt-f;,) 

-\^Jt2i [{d-^nrdAv) (a-V2t*"M;,) - {d-'flf^dAl,) (a-V2t*"5A2')] 

-iJtAid-'nf'dAy) (a-V2r"M;0 (3.34) 

4 Conformal symmetry on the light ray 

The conformal transformations act on the light ray positions z as 

, az + h 

z ^ z = 

cz + a 

This action is generated by 

^(0),- ^ gm ^ ^d, = -z'^d. (4.1) 

The generators obey the sl{2) algebra, 

[S°, S^] = ±S^, [S+, S-] = 25° (4.2) 

We need other representations i on functions of z generated by 

^W-^^(o)-^ Sii)^^z-'S^'^'/, S^'^+^z-''S^'^+z''. (4.3) 

The representation space V^^^ is spanned by z"^, and the function 1 represents 
the lowest weight state. 

The tensor product representation, V^^^^ <^V^^^\ is represented by polynomial 
functions of two variables '0(^1; ^2) and is generated by 8^2 — Si^^"' + 5'*^^^^". It 
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decomposes into irreducible representations, the lowest weight states of which, 

^- ^ = 0, 5°2^ = AV', are i)^^'^ = {zi - z^Y , with A = £i + 4 + n. 

The two-parton Hamiltonians obtained from QCD can be considered as acting 
on the tesor product representation spaces 

The operator is symmetric if 

(5f + 5f H = H (5^'^" + St^'^"") . (4.4) 
Writing the operator in integral form, 

Hi/j{zi,Z2) = j dZvdZ2i J{zi,Z2\zi',Z2')i>{zi',Z2'), (4.5) 

this results in the following condition on the kernel, 

(^f + Si''^" + + 4'"'"^") A^i, ^2^1', Z2') = 0. (4.6) 

Here wc assume that the integration is over closed contours in order to have no 
boundary terms in the interation by parts. The condition on the kernel is solved 
by the simple expression 

<''^i,£2A'A'v^i> ^2Fi', 2;2'j — 2;i2 ^i'2' ^ii' ^22' ^12' ^i'2 > l^-'J 

where the exponents are given in terms of the conformal weights ii and two 
paramenters d, M as 

1 1 
ai2^ d+1- - ai>2' = d-l + - 

ai2' = -^(^1 + ii' -k- k') + M, av2 = +^((^1 + h' -h- h') + M 
Oil' = -(-^1 - £i' + 1) - (i - M, 022' = (£1 - ^1' + 1) - - M 

a = 4 -4, +4-4' (4.8) 

The general solution is a superposition of the above expression with varying 
paramenters M and d. In order to close the contour one should choose a double- 
loop enclosing two branch points as in Fig. 2 (Pochhammer contour). For our 
applications we need integral operators where the integration is along the real 
axis in z'^ and z^- 

The kernel for closed contours transforms into the kernel for integration along 
the light ray by contraction of the contours to the branch cuts and rewriting the 
integral over finite intervals in a representation. 

In the regular case, if the exponents are larger than -1, we obtain the kernel 
in the form 

Z12 / dai / da2 S(zn' - ol\Zyi)^{z22' + a22;i2) 

JQ Jo 
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Figure 3: Double- loop contour. 



Q,«ii'Q,a22' (1 _ ^^p,, (1 _ ^^y,,, {1-ai- a2p'^' (4.9) 

Singular cases arise if some powers approach negative integers, e.g. if a = — 1 + £ 
then the kernel is obtained from the one above by by the formal substitution 

^ 5{a) + ej^. (4.10) 



^ \a 



Now it is not difficult to check that the kernels obtained in QCD have the con- 
formally symmetric form with the weights ^ = | for gluons and i = 1 for quarks. 

For example if ii = £2 = (-v = £2' = ^, M = 0,d = s we have in the integrand 
(4.9) 

(1 - ai - «2)^+^'-' / (i_«2)2^-i ^, (l-a,Y'-^\ 
^ ^ ■ d{ai)S{a2) + e {5{ai)- — p-^ + S{a2)- ^ ' 



(4.11) 

The second term results in the QCD kernel 

J?v~'\zuZ2\z,>,Z2') + (11' ^22'), 

•^11' = / da— — —[S{zu' - azu) - 6{zu')] S{z22')- 
Jo a 

encountered in the results (3.30 - 3.34). One checks also easily that the functions 
{zi — Z2)"' representing the lowest weight states are eigenfunctions of this kernel. 
The eigenvalues are 

Xn^ da^ ^ = V(l) -V'(2^ + ^) (4.12) 

Jo a 

They are proportional to the anomalous dimensions in the parallel helicity ex- 
change channel related to tranvcrsity parton distributions. 

The remaining kernels in the QCD operators correspond to the following 
particular cases of the conformal ones: 

99 99 (3.31), 



£1 = 4 = h' = 4' = 1 
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M — 0,d — —1 results in Jiu and M — 0,d — —2 results in J221. 

ff^ff (3.31), 

il=i2 = il' =i2' = l 

M = 0,d = —1 results in Jm and M = 0,d = —2 + e results in Jq\ 
ff^gg (3.34), 

^1 =4 = 1,^1' =4' = 1,(7 = -1 

M — 0,d = — I results in 2;J~2^ Jm and M — 0,d — + £ results in z^-^Jo^^- 
gg ^ //, (3.34) 

ii = h = iji' = i2' = i,cT = +1 

M = 0,d = — 1 results in 2:12^112 and M = 0,d = — | results in 2;i2'/22i- 
gf^gf (3.33), 

M = ^,d = results in Jm — | J211 and M = |, d = — | + £ results in 

7(9) , j(f) 
Jill -r J22' ■ 

gf ^ fg (3.33), 

^1 =^2' = 1,^2 = ^1' = 1 

M = 0, (i = — I results in J121 and M = 0,d = + e results in 



5 Integrable systems of interacting parton 

Consider operators acting on the tensor product (8) V2 V3 and in particular 
ones acting non-trivially on two of the three spaces (as indicated by subscripts) 
and obeying the Yang-Baxter equation (YBE) 

Ri2{u - v)Ris{u)R23{v) = R23{v)Ri3{u)Ri2{u-v). (5.1) 

R- operators obeying this relation can be used to construct integrable quantum 
systems. Consider N subsystems, the quantum states of which span the repre- 
sentation spaces V^^'\i — 1, in a periodic chain-hke configuration where the 
interaction betweeen the subsystems are due to exchanges along the chain with 
intermediate states out of the auxiliary space V^^°\ A complete set of commuting 
operators is obtained from the transfer matrix 

tM = Spo {Ri'r'^\u)Rg°''^\u)...Rit''-\u)) (5.2) 

by expansion in the spectral paramenter u. In the homogeneous case, ii — i,i — 
1, ...,N, these operators appear as sums of local operators (involving only a few 
neighbouring sites) for £0 — £. The first non-trivial operator in the expasion can 
be considered as Hamiltonian 

A.i^te{u)\^=o = E^v+i> H12 = P12 -^R^'''''\u)\^=o (5.3) 
Here P12 denotes the operator of permutation of subspaces. 
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The best known example of a R-operator is the one related to the Heisenberg 
XXX spin chain and is given by the 4x4 matrix 

= hx4 + {u + ija-^ ® a« (5.4) 

In this case the representation spaces correspond to spin ^ ( ^ = ~| in our 
notation) and are two dimensional. The tensor product with the unit matrix 
acting on the third subspace is not written cxplicitely here. 

The extension by sl{2) symmetry leads to another solution, called also Lax 
operator, appearing in YBE for one subspace of arbitrary representation i and 
two of the fundamental one. 

r[2'"'\u) = /2X2 ® I^'^ + {u+ ^)a" 4« ^{u+ + R2. (5.5) 

The goal is to construct the universal R-operator: R'^^'^'^'^\ which appears in 
YBE with ^3 = — I and ii , £2 arbitrary. This particular YBE implies the defining 
condition conditions 

(^1 + R2)R^'''''\u) = R^''^'''\u){Ri + R2): 
^{R^ - R2) + R,R2) R^'''''\u) - R^'^''^\u) (|(i?i - R2) + RiR2) (5.6) 

The first condition is equivalent to the condition of conformal symmetry which we 
solved before. The second fixes the freedom left in (3.8) for conformal operators. 

We look for the universal R operator in integral form. The conditions turn 
into differential equations for the kernel R{zi, 2^2 Z2')- 

(4'^) + Ri''^ + R^-''^ + R^'^') Rizu Z2\zr, Z2,) = 
^(i2f ) - )) + ^f^)^'^) + |(^S|-^^) - Ri)-'^^) + Ri'r'^^Rt''^) R = (|5.7) 

Here Rf^ is the non-trivial piece in the Lax operator matrix acing on Zi. 



R^ = _^wo l = (zj^'^ -^)^+^ 



; (^^ + ao) 



The decomposition of the 2x2 matrix structure into simple tensor products 
simplifies essentially the solution of (5.7). As an equivalent representation we may 
take the matrix element with repect of two vectors paramentrized by positions zi 
and Zr- 

< l\Ri\r >= {zi - l)Ri ( ] ) = 

~ZiiZirdi — £i{Zir — Zii) — Zj^j. Zi^~^ di Zj^^Zij^ ' (5-8) 
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We rewite the conditions (5.7) in terms of those generic matrix elements and use 
the fact that some derivative terms vanish for if zi or z.^. tends to some of the 
points Zi. In this way one obtains that the conditions (5.7) are equivalent to four 
first order equations. 

Vi R{zi, Z2\zi>, Z2') = 0. (5.9) 

Their consistency implies that the first order differential operators are similarity 
transformations of the simple derivatives, 

— Rd^R ^, 

and this transformation is just given by the wanted R operator kernel, 

1 ^2 1 Z\i , 2:2' ) — Zl^ Zy^ Z]^ Z]^ . (5.10) 

It is indeed a particular case of the conformal symmetric kernels with — ^2 = 
£2' , d — u, M — —u — 1, 

ai2 = u + 1- £i- £2, ai'2' = M - 1 + ^1 + 4 
ai2' = -li - 1 - £1 + £2, ai'2 = -li - 1 + ^1 - -^2 

Oil' = 022' =(7 = 0. (5-11) 

In the special case £1 = £2 we have 

jl{e,e) ^ I (^Q, / da25{zir - aiZi2)S{z22' + 0i2Zi2) 
Jo Jo 

(1 - ai - a2)"-'+''(l - ai)-"-^(l - aa)-""' 
We expand for small u using (4.10) 

S{Z2V)S{Z12') + 

u / da^^ 5(^12' - azu)d{z2v) + (11' ^ 22') + 0{u^) 

Jo [q;J + 

Comparing with the formula for the integrable Hamiltonian (5.3 ) we see that the 
terms with in the QCD operators (3.30 - 3.34) coincide with this Hamilto- 

nian for ^(/ = I and if = 1. 

It is remarkable that particular cases of parton interaction operators coincide 
with the integrable Hamiltonian of a generalization to (non-compact) representa- 
tions £ = 1 or I of the periodic homgeneous XXX spin chain. The multi-parton 
exchange contribution for parallel helicities of only gluons or only quarks is de- 
termined just by this operator. In the approximation of a large number of colours 
Nc closed chain configurations of gluons dominate and the corresponding con- 
tribution to the Bjorken asymptotics is calculated from the solution of the XXX 
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integrable system. In particular the lowest energy eigenvalues determine the 
powers of the large scale s of this contributions. 

Integrable systems corresponding to other helicity configurations of QCD par- 
ton interactions have not been constructed yet. There are some results related to 
open chains with three sites [18] . Besides of the terms related to the integrable 
Hamiltonian J^f,^^ + (1 2) there are other terms in the QCD results. However 
there is a close relation of the other operator terms appearing in QCD to the R 
operator. 

For the gluon case (3.31) Pi2-R''^'^^('u) for u = — 1 results in the kernel J112 
and for u = —2 in J22i- For the fermion case (3.32) Pi2R^^'^\u) results for u = —1 
and —2 is Jm and Jq'\ The kernels of the annihilation interactions (3.34) are 
obtained as PuR'^^'^^u) with u — —1,-2 for // — > gg and zi2Pi2R^^''^\u) 
with u — —1, —2 for gg f f. 

In case of mixed quark-gluon interactions Pi2R^^''^\u) with u = — | re- 
produces J211 and by crossing also Jm — |Ji2i and with m = — | reproduces 

j[y — J^y — . In this case there is one contribution left, the kernel + 
not reproducible from this R-operator kernel. In the supersymmetric extension 
this kernel appears as a component of the interaction Hamiltonian of XXX chain 
constructed from the superconformal sl{2\l) solution of the Yang-Baxter equation 
[21]. 
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